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Abstract 
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We propose a classification of critical behaviours of branched polymers for arbitrary 
topology. We show that in an appropriately defined double scaling limit the singular 
part of the partition function is universal. We calculate this partition function 
exactly in the generic case and perturbatively otherwise. In the discussion section we 
CN \ comment on the relation between branched polymer theory and Euclidean quantum 

O ! gravity. 
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1 Introduction 



The statistical mechanics of branched polymers (BP) is one of the simplest and 
most tractable models of random geometry. It is a subject of intrinsic interest 
and has already been studied by several authors A further motivation for 

developing these studies is provided by the recent suggestion || that important 
features of Euclidean quantum gravity can be inferred from those of the ensemble of 
branched polymers isomorphic to the ensemble of trees of baby universes connected 
by wormholes. 

The main body of this paper is devoted to the discussion of the critical properties 
of nakedf]] BP models, including the possibility of loop formation. Some results have 
been obtained earlier by other authors and are included here for completeness. We 
shall comment about quantum gravity in the discussion section. 

We formulate the problem as a minifield theory, defining a generating function 
W(j) by the familiar equation: 

e AW U ) = J d4>e M-^+vw+^] (1) 

However, here the integration variable <ft is just a rea l number. The interaction 
potential is assumed to have the form 

V{4>) = E(^A)/ ,k > 2 (2) 

k 

with positive couplings: pk > 0. Thus, strictly speaking, the integral in (D) does not 
exist. It is introduced in the first place to define a perturbation series. Obviously, one 
can associate diagrams with the terms of this series. The propagator corresponding 
to a link is A -1 , a source j is attached to each external link and the number of loops 
in the diagram is identical to the power of A -1 . The generating function of rooted 
diagrams (those with one marked external point) is given by dW/dj. 



2 The trees 

The saddle point = Z is found from the equation 

z = \- 1 [ J + Y,PkZ k - 1 ] (3) 

k 

When calculated in the saddle point approximation, the functions W = W and 
dW /dj generate the tree and the rooted tree diagrams respectively. Since 

Z = dWo/dj , (4) 

it is evident that Z is the partition function of rooted BP. And indeed eq. (Q) is 
generally used as the defining equation for rooted BP. We write here j instead of p\ 
and we have, for definiteness, incorporated P2 into A. But this is just a matter of 
conventions. 

3 That is without matter fields living on them. 
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Introduce the positive definite function 



W) = 4 + E^ fc ~ 2 (5) 

and rewrite (|3|) as 

A = F(Z) (6) 

The positivity of pk implies that F'(<p) can have at most one zero for <fi > and that 
F"((f)) > 0. Let r denote the radius of convergence^ of the series in (^]) and let <fi = <f) 
be the point in the interval (0,r] where F((p) takes its minimum value. We deduce 
from (|6|) that the parameter A cannot decrease below A c = F(<f)). Consequently, the 
partition function Z must have a singularity at A = A c . 

Write 5X = A — A c for later convenience and assume that in the neighbourhood 
of 5X = the singular part of the partition function behaves as 

Z s ~ ^A 1 ^ , (7) 

This is the conventional definition of the (geometrical) susceptibility exponent 7. 
Another interesting quantity is the two-point correlation function C(x), where x is 
the (integer) distance between two marked points. It has been shown in that at 
the tree level this function is 

C(x) ~ \\- x ^{ZF{Z)r (8) 

For large x 

x logC(x) = - const 5\ d H (9) 

where du is the Hausdorff dimension. 

Generalizing the considerations of ref. one can present the catalogue of pos- 
sible critical behaviours: 
The generic case. 

In the generic situation <j) < r. One then has F'(4>) = and 

F{Z) = F{$) + V(0)(^-0) 2 + ..., (10) 
which after inversion yields^ 

[2/F"(0)]3 5A5 (11) 

Hence 7 = ~. One also finds dn = 2. 

The alternative to the generic case occurs when = r < 00. It follows from 
the positivity of pk that (ft = r is a singular point of F((j)). Since F(<f>) decreases 
monotonically towards F(r) > one expects the singularity to be a branch point in 
all cases of physical interest^. Hence (modulo logs) 

F ( Z ) = E ^-^F [k \r)(r - Zf + c(r - Zf~ x + ... , n + I < (3 < n + 2 (12) 
k - 



4 Possibly infinite, but not infinitesimally small. 

5 We have written a minus sign in front of the 2nd term on the right-hand side because we 
consider the branch Z < <j> as the physical one. 

6 The interested reader can consult || and in particular the theorems by Leau, Le Roy and 
Lindclof. 
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Since F is concave, one must have (3 > 2. Furthermore, from the relation between 
F and V and from the positivity of all the derivatives of V one easily deduces that 
c(— l) n < 0. There are two possibilities: 
The semi-generic case. 

It occurs when = r and F'(r) < 0. Inverting ([TJ) one finds 

5\P~ l (13) 



[-F'{r)\ 



and therefore 7 = 2-/5. In this case du = oo. 
The marginal case. 



It occurs when <p = r and F'(r) = 0. For 2 < /3 < 3, inverting eq. (|12|) one finds 



Z S ^(-)^ (14) 
c 

so that j = (J3- 2)/ (f3 - 1). The Hausdorff dimension is d H = ((3 - 1)/(J3 - 2). 
When (3 > 3 the situation is analogous to the generic one: one has ([TTJ) with r in 
place of and 1 = \ while d H = 2. 

To summarize, in the non-generic cases there is a continuum of universality 
classes, characterized by the exponent (3 and by the vanishing (or non- vanishing) of 
F'(r). Each class corresponds to an infinity of different choices of the couplings pk- 



3 The loop expansion 

3.1 The generic case 

The loop expansion is obtained calculating corrections to the saddle-point approx- 
imation. This can be done either directly or by using Dyson- Schwinger equations. 
We employ here the second method, which is more elegant and enables one to write 
rapidly the BP equation, an analogue of the string equation one finds in the double 
scaling limit of 2d gravity |T0|. 



Set U(cb) = |0 2 — V(4>). The Dyson-Schwinger equation reads 

[^(A- 1 |r)-j]e Awr W = (15) 

which can be rewritten as 

U\Z + K- l ^.)-l=3 (16) 

Here Z = dW/ dj denotes the full partition function, with loop corrections included. 
Let 0o be the point where U" vanishes: U"(<fto) = 0. Define 

A = 2[j - ^(0 O )]/C/ W (0o) (17) 

and the double scaling limit: 

A^0,A^oo,t = -i U"'((bo) A At = const (18) 
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Both A and t are positive definite. Set Z = <fio — A 2 x(t) in ( P^D and expand. In 
the limit (|18D one gets the BP equation promised at the beginning of this section: 

l = X 2 (*)+^x(t) +*'(*) (19) 
From the Riccati equation (|T9| ) one can obtain x as a universal power series in t" 1 : 

X(t) = 1 - — - -At + ... (20) 
AK 1 6t 72t 2 K ' 

We have assumed that the first term is +1. This is the physical choice (cf. the 
footnote preceding eq. (|TT1)), which also guarantees that the terms corresponding 
to higher topologies give a positive contribution to the partition function. The 
behaviour of the coefficient \ n of t~ n , for n ^> 1, can be easily estimated using (|19l): 

F(n) , . 

Xn ~ const (21) 

The constant above can be found numerically and is approximately equal to —0.32. 
As one could expect, the series ( |20"D is not Borel summable. The perturbative series 
does not determine %{£) uniquely. The first singularity of the Borel transform occurs 
when its argument equals 2. Thus the leading non-perturbative contribution to 
is proportional to e~ 2t . 

Let f(t) be the primitive of f'{t) = x{t)- Setting $(t) = e-^ and changing 
the independent variable t —> z = (§£)» one gets from ( ]TP| ) the Airy equation 

$" = z§ (22) 

Hence 

$(2) = a Ai{z) + b Bi(z) (23) 

and 

vf^ = ( . , - - , • 

,2 7 $(z) ' V2 

The function x{t) depends on a single, but arbitrary parameter a/b, which measures 
the strength of the nonperturbative contribution to the solutionf]. In this respect 
the situation resembles that of string theory. 

It remains to find the relation between A and SX. The position of the point 
= 0o depends on A, but not on j. One easily finds, however, that as A approaches 
Ao, its j-dependent critical value, one has 0o — > 4> + 8X/V'"(4>), which is also j- 
dependent. One further finds from the definition of A that 

2 ° 5\ [1 + 0(5X 2 )} (25) 



1 T ,/ s „ 2 

3 



, (t) = f5 t V'|M, z = f3 t) (24) 



V"'((t>) 

The behaviour of x{t) at large t and the above equation imply the following depen- 
dence of the susceptibility exponent on the number of loops L: 

^=\ + \L (26) 

a result obtained by a different method in ||. 

7 In the limit b — > one obtains the unphysical solution with x(°°) = ~ !• 
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3.2 The strong coupling regime 

Set V(4>) — 9V(4>) and assume that g 1. This defines the strong coupling regime. 
It is obvious from @ that in this regime ~ \j 3 1 9 an d ^"'(0) ~ 9- Hence, for 
large enough g one necessarily has < r and one is in the generic situation. One 
easily finds 

.3 

t ~ ^ A 6\% (27) 
9* 

Therefore, in the double scaling limit t remains a small parameter. Using ( p4j) one 
gets 

. i 

Z si = -VK X (t) ~- J 4V5\ (28) 

9 s 

The free parameter a/b is hidden in the coefficient multiplying the right-hand side 
of fl28P . The susceptibility exponent 7 = § and the number of loops is zero0. 

3.3 The non-generic cases 

When the minimum value of A is found at the boundary of the convergence interval 
of the potential, the equation U"(4>o) = has no solution with 0o < r - in this case 
we combine the Dyson- Schwinger equation ([16]) with the expansion ([12]). We then 
encounter formal expressions of the type [h(x) + d/dx}^" 1 , which must be given a 
precise meaning. Write 

Mx) + d/dxf' 1 = — / ds s -P e -»lHx)+d/dx) (29) 

r(i — p) 3o 

Define G by 

e -s[h(x)+d/dx] = e G e -sd/dx / 3 qn 
It is easy to check that e G satisfies the differential equation 

4-e G = -e G e sd/dx h(x)e- sd/dx = -e G h(x - s) (31) 

OS 

with the boundary condition G — 1 at s — 0. Hence G is a c-number function: 

G — H(x — s) — H(x) (32) 
where H(x) is the primitive of h(x): H'(x) = h(x). Using this result we obtain 

[h(x) + d/dxf-' ■ 1 = ± ^^ h^RM (33) 
where Rk{h) is defined by the equation 



^ ( ^ ) = ^^ [e " (X " S)e " W] ^ ( 34 ) 



3 The average number of loops equals 9 In Z/dlnN calculated at fixed SX. 
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and depends on the derivatives h'(x), h)- n ~ 1 ^ (x). Using Leibniz formula to calculate 
the derivative in fl34|) and replacing d/ds by d/dx in the appropriate place one can 
get rid of s. After some algebra one obtains the following recursion formula 

Mh) = \\h\x)R k „ 2 {h) - R' k _M] (35) 



with R = 1 and R\ — 0. As one might expect the series ( p3|) would be truncated 
at k — (5 if (3 were an integer (which it is not!). We shall now use the above results 
to study the non-generic scaling behaviour. 
The semi-generic case 
We define 

A = U'(r) — j = rSX (36) 
and define the scaling limit as follows: 

AA 2 

A -»• , A -> oo , t= — - = const (37) 

—rF\r) 

Rewrite the Dyson- Schwinger eq. (16) in terms of F, set Z = r — f and use the 
expansion (|T2| ) to get 

= A + rF'(r)f - rc(f + -^Y' 1 ■ 1 + - (38) 

where the dots represent terms irrelevant in the scaling limit. 
We have f = f a — Z s , where 

/a = - (39) 

—rb \r) 

is analytic in A. In the case under consideration, Z§ ~ A' 3 ^ 1 is subleading, compared 
to fa, because (3 > 2. But we are precisely interested in this subleading term. 
In the scaling limit one gets from (f 



^zj^ + Js)"- 1 (40) 

We can now use (|33|), slightly modified: because of the factor A -1 in front of ^ one 
has R k (f) = A~ k R k (Af) instead of R k (f): 

zs—^^^f^-nw) (4i) 

One finds by inspection that in the scaling limit only the terms with even k con- 
tribute to the right-hand side. The leading contribution comes from the first term 
in 

R2m = C 2m (fV + - (42) 

We find from the recursion equation (|35D that 



Using (fjgj), (gg) and flEJ) we finally getfl 

From the above equation one reads the susceptibility exponent 

lL = 2 - p + 2L (45) 

The marginal case 

One again has the expansion (|5&D , but with F'(r) set to zero. For /3 > 3 the 
problem reduces to the generic one: it is sufficient to replace in all formulae cf) by r. 
When 2 < /? < 3 the scaling limit is defined as follows: 

A -> , A -> oo , t = (rc)~A(A) " 1 = const (46) 
In this limit the problem reduces to the solution of the equation 

- = (/ + tI^- 1 • i (47) 

rc AaA 

with / = — Z s = (A/rc)~x(t)- We have not succeeded in calculating in a closed 
form the asymptotic expansion of x{t)- What we can offer is a systematic recursive 
scheme, more efficient than eq. (|33|) . Before entering into further algebra and 
anticipating on the result to be obtained, let us mention that the expansion is in 
inverse powers of t, which implies that the susceptibility exponent is 

7i = ^T + ^T (48) 

Introduce an auxiliary variable x = t~ and rewrite (17) as 



x 



r(i-p) Jo 



d s s -P e G(s,*) (49) 



where G is given by (32) and h(x) = t^x(t)- Define u = t 1 and perform the 
substitution 



i 



s = aw 3 (50) 
x = u^ 1 (51) 

After simple algebra one obtains a differential equation with analytic coefficients 
r . 1 an, dG 2 dG .„ 1. . . 

l(1 -^ ) + T 1 *-" ^ = - {1 -p )x (52) 



9 Remember that for n + 1 < (3 < n + 2 one has c(— l) n < 0. This insures that the terms with 
2m > — 1 are positive. 
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We write G and x as power series in u and insert these expansions into eq. fl52|) to 
get a hierarchy of differential equations for the coefficients G n (a). The solution to 
these equations is Go(a) = —a (since xo = 1) and 



n+i r(f4 + fc) 

" -I 1 g 

-r(A; + f)r(^ + i; 



G n (<r) = ~ E w,, , ,wn-fc , ^ ° k Xn+i-k ,n>l (53) 



3 

Notice that G n depends on Xk with k < n only. Finally, eq. ( |4"9"D becomes 

Equating to zero the coefficient of t _ra on the right-hand side yields Xn in terms of 
{ Xk ■ k < n }. 



4 Discussion 

Let us summarize what has been achieved in this work. We have generalized the 
discussion of ref. [H], showing that BP models fall into one of three categories. 
Furthermore, we have extended the discussion to arbitrary topology. For each of 
the categories in question we have defined an appropriate double scaling limit and 
written the singular part of the partition function as a universal asymptotic series, 
each term of the series corresponding to a given topology. In the generic case, we 
have derived a BP equation, valid also in the non-perturbative regime. In short, we 
have done in the context of BP models what has been earlier achieved for matrix 
models. 

As in the latter case, one can relax the constraint > and consider multi-criti- 
cal model, where there exists a point <\> = 0o such that U"((po) = ... = [/[ m_1 ] (0 O ) = 0. 
One obtains an equation analogous to (0). However, now one has the integer m 
instead of f3 — 1, so that the equation is a genuine differential equation of order 
m — 1. Using the techniques developed in the last section one defines a universal 
expansion of the partition function 

oo 

Z s ~ -A£(l - Xnt~ n ) , t ~ NA^ (55) 
=i 



However, the sign of the coefficients Xn is n ot positive definite, as expected. 

Our motivation for studying this problem originated from our involvement in 
the study of simplicial gravity. Therefore, we would like to end this paper with 
remarks concerning the hypothetical relevance of the study of branched polymers 
for quantum gravity. 

One of the striking results obtained from computer simulations of random ge- 
ometries is that a random manifold does not stay smooth during the simulation. It 
develops a tree structure: the nodes are the baby universes and the links are the 
bottlenecks (wormholes) that connect them. We shall call this tree the skeleton tree 
of the manifold. For example, in 4d one observes a transition between two phases. 
In one of them the skeleton trees resemble generic BP. In the other the manifolds 
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are crumpled: there is one big mother universe and a large number of small babies 
attached to it. Similar phase transitions are encountered in lower dimensions. 

It has been suggested in ref. || that these phase transitions reflect the dynamics 
of the skeleton trees. Indeed, in BP theory the position of the minimum of F{(f)) 
changes when one moves in the coupling space and a phase transition occurs when 
this minimum hits the boundary of the support of this function: the generic BP turn 
into crumpled structures i.e. with large or infinite Hausdorff dimension. In sect. 2 
the Hausdorff dimension has been calculated at the tree level, but the result seems 
to hold for any fixed topology. 

The hypothetical relation between the physics of 4d manifolds and that of the 
associated skeletons, if confirmed by further studies, may turn out to be a fruitful 
idea. The present numerical studies of 4d random manifolds are limited to simplest 
topologies. Thus the models under study are non-unitary: the wormholes that have 
been emitted cannot be reabsorbed. These simulations, although limited to skeleton 
trees without loops, are however sufficient for the determination of the couplings at 
the nodes (the s). A glimpse at the unitarized theory is provided by an extension 
to trees with loops, along the lines of this paper. 

Notice also, that a phase transition analogous to that observed in 4d occurs in 2d, 
but only for large enough c (or at least for c > 1). But a sensible 2d gravity theory 
exists for c < 1 onlyQ Thus one is led to speculate that a sensible gravity theory in 
4d will only be obtained when one does in 4d something analogous to the reduction 
of the central charge in 2d. One can develop a heuristic argument analogous to that 
proposed by Cates [0 to explain the c = 1 barrier in 2d. Consider the continuum 



theory with Einstein-Hilbert action and write the metric in the form g a b = e 2a g a b. 
Assume that the BP phase in 4d is dominated by the dynamics of the conformal 
factor and use the effective action calculated in [|I4J to estimate the free energy of a 
diluted gas of "spikes" : e 2<J = 1 + p 2 /[{x — x ) 2 + a 2 ], with a<p< const x a. As 
discussed in (|Tj| the parameter a can be taken arbitrarily small, while keeping the 
invariant cut-off fixed. The free energy is 



_ r 1411 + N s + 11N F + 62N V -28 , 1 

F ~ 4 In - , 

L 360 J a ' 



(56) 



where 1411 is the contribution of transverse gravitons, Ns,f,v is the number of 
scalar, fermion and vector fields respectively and —28 is the quantum a and ghost 



contribution. Notice that, as emphasized in |14| , the sign of matter and ghost 
contributions is in 4d opposite to that found in 2d. For pure gravity the coefficient 
in front of In - is negative and one expects a condensation of "spikes". However the 
addition of matter fields can easily change the sign of this coefficient and stabilize 
the theory. 
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10 This is at least true for generic random surface models. Exceptions to this rule have been 
claimed [ p~T| . 
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